A Cluster of Compact Radio Sources in W40 
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ABSTRACT 

We present deep 3.6 cm radio continuum observations of the HII region W40 obtained using 
the Very Large Array in its A and B configurations. We detect a total of 20 compact radio sources 
in a region of 4' x 4', with 11 of them concentrated in a band with 30" of extent. We also present 
JHK photometry of the W40 cluster taken with the QUIRC instrument on the University of 
Hawaii 2.2 meter telescope. These data reveal that 15 of the 20 VLA sources have infrared 
counterparts, and 10 show radio variability with periods less than 20 days. Based on these 
combined radio and IR data, we propose that 8 of the radio sources are candidate ultracompact 
HII regions, 7 are likely to be young stellar objects, and 2 may be shocked interstellar gas. 

Key words: HII regions - infrared: stars - open clusters and associations: individual (W40) 
- radio continuum: stars ~ stars: formation 



1. Introduction 

The W40 HII region (denoted as S64 in the 
Sharpless catalog) is the central object in a rich 
cloud complex that lies some 3.5° above the 
Galactic pla ne at J2000 coordinatesJ^SF^^Sf, 
-02°05'4" (IWesterhoud [1958 : Sharplesd Il959l) . 
W40 is a blister HII region with a diameter of 
'--'6'. It is bordered by a molecular cloud that 
subtends an angle of ~ 1°. The dense core of this 
molecular cloud is located at galactic coordinates 
28.77 +3.70, and is labe led TGU 2 79-P7 in the 
all-sky cloud atlas of iDobashi et al.l (|2005l ). 



W40 has been well-studied at radio and mil- 
limeter wavelengths, leading to a firm understand- 
i ng of the region's gaseous structures. ICrutcher fc Chu 
(jl982l ) compiled a comprehensive review of the 
available molecular line and radio continuum data 
to develop a detailed kinematic picture of the 
HII region and its interaction with the molec- 
ular cloud. The warm CII interface between 
these two principal components of the clou d com- 
plex has been extensively studied as well ( Vailed 
19871 : IVallee fc MacLeodlll99ll : IVallee et al.lll992[) . 



Additionally, Zeeman measurements of OH ab- 
sorption lines have revealed a magnetic field of 



B=-14.0 ± 2.6/iG (ICrutcheret al.lll987l) . 

A dense stellar cluster within W40 is apparent 
in th e IR maps of the 2M ASS and DENIS sur- 
veys (JRevle fc RobinI 120021) . This cluster is dom- 
inated by three central OB stars that constitute 
the primary excitat i on sources of the HII region 
(jZeilik fc Ladall978HSmith et al.lll985l ). The stars 
in this cluster appear to be heavily obscured along 
our line of sight , with Ay ^ 10 mag nitudes of vi- 
sual extinction (jShuping et al.lll999t ). W40 IRS2a 
has been suggested as the dominan t source of ion- 
izing radiation in this region ( Smith et alJll985) . 
and millimeter observations have shown evidence 
for circu mstellar dust shells aro und IRSla.. 2a, a nd 
3a (Smith et al1[T985; Vallee fc MacLeodlll994| ). 



From the collective results of these studies we 
can see that the W40 HII region has broken 
through the surrounding molecular cloud in a di- 
rection that lies at an angle to our line of sight. 
Thus, the stellar cluster which powers the HII re- 
gion is mostly hidden in the optical and near-IR, 
but scattered Ha emission marks the blistering 
edge of the HII region. The dense core of the 
molecular cloud is offset from the HII region and 
its embedded cluster by '^ 2', and a thin CII region 
delineates the boundary between the cool molecu- 
lar gas and the hot ionized region. 

The distance to W40 is not yet well constrained. 
Several measurements of atomic and molecular 
lines along the line of sight towards W40 have 
been used to estimate the distance to W40 — as- 
suming that the cloud complex i s in a circular or- 



bit a r ound the gala c tic center (IReifensteiii et a.1. 
19701 : iDownesI [l970l : iRadhakrishnan et al.lll9m 



Taken together, these measurements produce only 
a weak constraint of 300-900 pc. Using assump- 
tions about the spectral types of several of the 
stars in the W40 cluster, distance modulus es- 
timates can only limit t he distance to a wide 
range from 400 to 70 PC (ICrutcher fc ChiJll982l 



Smith et al.l Il985l) . iKolesnik fc lurevichl (figj 



calculated a distance of 600 pc using a novel dis- 
tance determination technique involving OH line 
measurements. Hereinafter we will adopt a dis- 
tance of 600 pc, which puts the W40 complex at 
a height of 37 pc above the galactic plane. 

More recently, the properties of the W40 
region have been s u mmar ized in a review by 
Rodnev fc ReipurthI (|2008[ ). Also, iKuhn et al ' 



and detected about 200 young stellar members of 
the cluster. 

2. Observations 

We have used the Very Large Array (VLA) 
of the NRAO in its B and A array configura- 
tions to observe W40 at 3.6 cm, under the NRAO 
program code AR551. The B configuration ob- 
servations were made on 2003 November 3 while 
the A configuration observations were made on 
2004 September 18. Our phase center was at 
a(J2000) = 18''31"23!7; (5(J2000) = -02°05'29". 
The amplitude calibrator was 1331 + 305, with an 
adopted flux density of 5.18 Jy, and the phase 
calibrator was 1804 + 010, with flux densities of 
0.76 ± 0.01 and 0.82 ± 0.01 Jy for the first and 
second epochs of observation, respectively. The 
half power contour of the synthesized beams were 
0;'88 X 0"76; +17° (for configuration B data) and 
0:'26 X 0"22;+16° (for configuration A data), for 
images made with the ROBUST parameter of the 
task IMAGR set to (Briggs 1995). The data were 
analyzed in the standard manner using the NRAO 
AIPS package. The data for the two epochs were 
self-calibrated in phase separately and then con- 
catenated to obtain a final image with the high- 
est sensitivity possible. The half power contour of 
the synthesized beam for this concatenated data 
was 0''36 X 0''30; +28° for images made with the 
ROBUST parameter of the task IMAGR set to 
(JBriggs 1995) . 

A total of 20 sources were det ected in a field 
of 4' X 4' in this image. Following iFomalont et al. 
()2002l ) , we estimate that in a field of 4' x 4' the 



a priori number of expected 3.6 cm sources above 
~0.2 mJy is only ~0.3. We then conclude that 
perhaps one out of the 20 sources could be a back- 
ground object, but that we are justified in assum- 
ing that all the detected sources are associated 
with W40. 

Infrared photometry was carrie d out with the 
QUic k Infrared Camera (QUIRC, iHodapp et al 



19961 ) on the University of Hawaii 2.2 meter tele- 
scope from 2003 June 17 to 19. Data were ob- 
tained in the J, H, and K' filters of the Mauna 



(J2010I ) have performed a deep X-ray study of W40, 



Kea O bservatories (MKO) filter set (jTokunaga fc Vacca 
20051 ). These IR images were reduced and com- 
bined using the Gemini Observatory's QUIRC 
package in IRAF. Our data were collected un- 



der non-photometric conditions, and errors in the 
telescope guiding combined with persistent focus 
problems made it impossible to pursue photome- 
try by fitting the point spread function, as is pre- 
ferred for crowded fields. Thus, aperture photom- 
etry was performed using the NOAO DIGIPHOT 
package. The four brightest IR stars were satu- 
rated even in our shortest (1 sec) exposures, so 
without PSF fitting we were only able to obtain 
upper limits on the magnitudes, which agree well 
with previous work. 

We find over 1000 IR sources with detections in 
at least one of the three bands, of which 271 have 
reliable photometric measurements in all three 
bands. A subset of 118 sources that have good 
JHK' photometry in our QUIRC data and good 
JHK photometry in the 2MASS catalog were used 
as calibration stars to translate the entire dataset 
to the standard MKO system. This calibration 
step is the dominant source of error in our IR pho- 
tometry, because it relies on transformation equa- 
tions which convert our measured K' photometry 
into the K band and take the 2MASS JHK mag- 
nitudes into the MKO JHK system (Conn elley, 
private communication, iLeggett et al.l 120061) . In 
Figure [T] we demonstrate the scatter in our cal- 
ibration step by plotting the J magnitudes we 
have measured and calibrated to the MKO system 
against the J magnitude reported by 2MASS, also 
translated to the MKO system. We estimate that 
this calibration procedure introduces about 0.1 - 
0.2 magnitudes of error into our final photometry 
values. A table of calibrated JHK photometry 
for all 276 well-measured IR sources is available in 
the online edition of the Astronomical Journal. 

Of the 20 radio sources detected in our VLA 
data, 15 have counterparts in our IR data with 
positions matched to within 1". In Table [2] we 
list the positions and flux densities of these 20 ra- 
dio sources, averaged over our two observations. 
We also note in column 5 if they were found to 
be radio variable or not, by comparing the ob- 
servations of the two epochs and restricting the 
(u,v) coverage of the images for comparison to the 
range of 16 to 300 kA, a coverage in common for 
the A and B configurations. The number given 
in parentheses for the variable sources is the ratio 
between the largest and smallest flux density ob- 
served. Columns 6,7, and 8 provide the measured 
J, H, and K magnitudes of the IR counterparts, if 
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Fig. 1. — Detected IR sources that appear in the 
2MASS catalog, showing our measured J magni- 
tude plotted against the reported 2MASS J mag- 
nitude. Both axes have been translated to the 
MKO filter system. 



any. As noted above, the four brightest IR coun- 
terparts (VLA-1,7,13, and 15) were saturated in 
our data. For these sources we have translated 
the published magnitudes in the 2MASS catalog 
to the MKO system using the same transformation 
equations employed for our calibration procedure. 
Typical internal errors on our measured IR values 
are <0.1 magnitudes. 

3. Overall Characteristics of the Sample of 
Radio Sources 

A composite JHK' image of the entire W40 re- 
gion is shown in Figure [2l This field encompasses 
a region of ~ 7' centered around the stellar cluster 
and the middle of the HII region. The two boxes 
in this figure indicate regions shown with greater 
detail in Figures [3] and 21 Figure [3^ displays a ra- 
dio contour map containing the entire radio cluster 
detected at 3.6 cm. Out of the 20 detected radio 
sources, 11 of them — including the 2 brightest 
— are concentrated in a band with 30" of extent, 
as shown in Figure |4^. Figures |3)d and |4l3 show 
the same regions in the K' filter, revealing that 



Table 1 
Parameters of the 3.6 cm VLA Sources 



W40-VLA 


R.A. 


Dec 


Flux 


Time 


J 


H 


K 


W40-IRS 


Number 


0^2000 


(52000 


(mJy) 


Variable? 




(MKO) 




Designation* 


1 


18 31 14.822 


-02 03 49.99 


0.92 


N 


8.22 


7.507 


6.98 t 


_ 


2 


18 31 15.279 


-02 04 15.19 


0.90 


N 


- 


- 


- 


- 


3 


18 31 22.328 


-02 06 19.60 


0.47 


N 


14.27 


12.17 


11.10 


- 


4 


18 31 23.244 


-02 06 18.04 


0.19 


Y(>2.0) 


14.27 


12.72 


11.72 


- 


5 


18 31 23.622 


-02 05 35.80 


3.97 


Y(5.3) 


11.92 


10.33 


9.70 


- 


6 


18 31 23.687 


-02 05 27.90 


0.66 


Y(>2.0) 


- 


- 


- 




7 


18 31 23.985 


-02 05 29.53 


2.40 


Y(3.2) 


8.64 


7.21 


6.06 t 


IRS 2a 


8 


18 31 26.021 


-02 05 17.02 


3.27 


Y(8.3) 


10.74 


9.09 


8.15 


IRS Ic 


9 


18 31 27.312 


-02 05 04.44 


0.99 


Y(1.6) 


13.51 


12.03 


11.17 


- 


10 


18 31 27.459 


-02 05 11.96 


0.82 


Y(>5) 


11.41 


9.68 


8.79 


- 


11 


18 31 27.585 


-02 05 17.68 


0.95 


Y(6.4) 


13.65 


12.18 


11.20 


- 


12 


18 31 27.641 


-02 05 13.46 


1.64 


N 


13.50 


11.68 


10.47 


- 


13 


18 31 27.657 


-02 05 9.68 


0.86 


Y(>5) 


10.57 


9.36 


8.59 t 


IRS Id 


14 


18 31 27.681 


-02 05 19.66 


5.78 


N 


11.49 


9.91 


8.46 


- 


15 


18 31 27.812 


-02 05 21.86 


1.71 


N 


7.60 


6.50 


5.53 t 


IRS la 


16 


18 31 28.016 


-02 05 17.90 


0.94 


N 


10.81 


9.44 


8.56 


- 


17 


18 31 28.646 


-02 05 29.04 


0.48 


N 


- 


- 


- 


- 


18 


18 31 28.654 


-02 05 29.74 


11.1 


N 


11.23 


8.93 


7.27 


IRS lb 


19 


18 31 28.683 


-02 05 22.25 


0.20 


N 


- 


- 


- 


- 


20 


18 31 30.213 


-02 07 18.03 


2.73 


Y(^17) 


- 


- 


- 


- 



*From lSmitli et al.l (|1985[ ). 

^Source was saturated in shortest exposure. Magnitudes quoted here are from 2MASS, translated to the 
MKO system. 
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Fig. 2. — Composite JHK' image of the W40 region. Two rectangles delineate the fields shown in Figures [3] and |31 



most of the radio sources can be associated with 
IR counterparts. 

Of the 20 sources detected, 10 were found to be 
time variable from one observation to the other 
(see column 5 in Table 1). This percentage (50%) 
of time variable sources is ve ry similar to t h e valu e 
of 47% found in Orion by IZapata et aP (|2004h . 
However, this percentage is large r than the value 
of 32% found in NGC 2024 by iRodn'guez et"aL 
(|2003l ). We searched unsuccessfully for circular 
polarization in the sample. We measured the an- 
gular size of the sources using the AIPS task IM- 
FIT, with a correction for bandwidth smearing. 
With the exception of sources VLA 8, VLA 13, 
and VLA 19, all sources were found to be unre- 
solved, Os < 0'.'2. 

4. Discussion 



There are several mechanisms that can produce 
compact centimeter radio sources in regions of star 
formation. In regions of low mass star formation, 
thermal jets and gyrosynchrotron emitters can be 
present. In regions of high mass star formation 
we also have strong ionizing radiation available 
and, in addition to the two mechanisms present 
in low mass star-forming regions, we can have ion- 
ized stellar winds, ultracompact HII regions, and 
radio proplyds. It is possible to distinguish be- 
tween these various possibilities with high angular 
resolution, multifrequency observations made with 
the required sensitivity. 

Let us first consider the subset of 15 sources 
that have IR counterparts. It is possible that these 
objects are ultracompact HII (UCHII) regions, 
centered around young massive stars. The radio 
emission can then be explained as optically thin 
free-free radiation, and the fa miliar photoioniza - 
tion equilibrium equations (e.g. IOsterbrocklll974l ) 
can be used to determine physical parameters of 
the system. However, UCHII regions are not ex- 
pected to exhibit short-period time variability, as 
is seen in 8 of the 15 radio/IR sources. The re- 
maining 7 sources can be considered as UCHII 
candidates. Among this group are the two bright- 
est radio sources, VLA- 18 and VLA- 14, as well as 
the brightest IR source, VLA-15 = IRSla. With 
an assumed electron temperature of 10** K and 
a distance of 600 pc we can estimate that a BO 
ZAMS star would be required to maintain the ion- 



ization in each of the se UCHII regions ([Panagia 
19731 : lThompsonlll984l ). 



As for the sources that show time variability, 
these are more likely to correspond to young stel- 
lar objects (YSOs) with active magnetospheres. 
The gyrosynchrotron radiation from such objects 
is expected to show substantial variation over 
timescales of a few days. One of the brightest 
radio/IR sources that may be placed in this cat- 
egory is VLA-7 = IRS 2a. The SED of this ob- 
ject has been examined from near-IR to millime- 
ter wavelengths, and it indicates that the star is 
surrounded by a pproximately 0.1 Mp) of circum- 
stellar material (jVallee fc MacLeod! [l993 ). This 
finding is consistent with our tentative classifica- 
tion of VLA-7 as a YSO, since many stars remain 
shrouded in their natal envelopes all the way to 
the main sequence. It is important to note that 
this classification is at odds with the results of 



Smith et al.l (|1985I ) and lVallee fc MacLeod! (|l991f ). 



who label this source as a deeply embedded OB 
star that is powering the ionization of the W40 
HII region. Based on the data presented here, the 
radio sources VLA-14, 15, and 18 may be equally 
valid candidates for supporting the HII region. 

Of the five radio sources that do not have IR 
counterparts, VLA-17 and 19 are of particular in- 
terest. VLA-19 is offset by 7'.'5 (0.02 pc) to the 
north of VLA-18, and VLA-17 is only 0'.'7 (0.002 
pc) away to the northwest. These radio sources 
are both relatively weak, do not vary with time 
on short periods, and have mildly non-circular ra- 
dio contours. All of this evidence suggests that 
these two radio sources may correspond to shock 
fronts from a thermal jet interacting with the am- 
bient interstellar medium. Based on its position 
and proximity, VLA-18 — IRSlb could well be the 
star which powers this radio jet. 

In order to unambiguously determine the na- 
ture of these compact radio sources, further infor- 
mation is required. A good indicator of UCHII 
regions is the spectral index at radio frequen- 
cies. Free-free emission from an UCHII region 
should have a very flat radio spectrum, so a spec- 
tral index more negative than -0.1 is a clear in- 
dication that the obje ct is not an UCHII region 
(JRodriguez et al.!!l993! ). This measurement could 
be made using the EVLA to observe W40 with 
similar angular resolution at 6 cm, for example. 
Another important observation to further our un- 
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Fig. 3. — (a) VLA continuum image of the compact radio cluster in W40 at 3.6 cm. The half power contour 
of the beam is shown in the bottom left corner. The contours shown are 5, 6, 8, 10, 20, 40, 80, and 160 times 
21 /iJy beam~^, the rms noise of the image, (b) Composite IR map of the same region, showing that most 
of the bright radio sources have IR counterparts. 
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Fig. 4. — (a) Radio contour map as in Fig. [3^. This close-up view of the central region shows 11 of the 20 
radio sources, (b) Composite IR map covering the same region shown in (a). 



derstanding of the W40 cluster will be to gather X- 
ray data to determine the high-energy character- 
istics of the dominant radio/IR sources. Gyrosyn- 
chrotron emitters should produce detectable X-ray 
emission, with measurable time-variability signa- 
tures. A combination of high angular resolution 
radio observations at longer wavelengths plus X- 
ray measurements with XMM or Chandra would 
enable us to definitively determine the astrophysi- 
cal nature of these sources. A first Ch andra study 



of W 40 is currently in preparation (jKuhn et al 
20in[ ). 



5. Conclusions 

We have carried out the first detailed, high res- 
olution study of the center of the W40 region at 
3.6 cm, supplemented by deep IR imaging of the 
entire W40 complex. Previous radio observations 
of W40 were not able to distinguish individual 
sources, but our data have revealed a cluster of 
some 15-20 radio sources. Many of these sources 
are clustered together in a narrow band at the cen- 
ter of the W40 cluster. A substantial fraction of 
these radio objects are also present in our JHK' 
images, and a number of them show significant 
variability over timescales of a few days to weeks. 
The combination of this high resolution radio data 
and our IR maps allows us to speculate on the na- 
ture of these sources, possibly distinguishing be- 
tween YSOs, UCHII regions, and shocked inter- 
stellar gas. 

W40 is a rich, nearby star forming region that 
has not benefited from the same intense scrutiny 
afforded to regions such as Orion, Ophiuchus, and 
others. This is largely due to the veil of gas and 
dust which hides most of W40's compact objects 
from view. Radio and IR observations such as 
those described here are necessary to penetrate 
the obscuring shell to investigate the interior. X- 
ray observations and further high resolution radio 
measurements at longer wavelengths will be very 
beneficial in unmasking the W40 cluster. 

This material is based upon work supported by 
the National Aeronautics and Space Administra- 
tion through the NASA Astrobiology Institute un- 
der Cooperative Agreement No. NNA08DA77A 
issued through the Office of Space Science. This 
research has made use of the SIMBAD database. 



operated at CDS, Strasbourg, France, as well as 
NASA's Astrophysics Data System Bibliographic 
Services, and the NASA/ IPAC Infrared Science 
Archive, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronau- 
tics and Space Administration. The infrared data 
presented here were obtained with the UH 2.2m 
telescope on Mauna Kea, operated by the Univer- 
sity of Hawai'i. The authors wish to recognize and 
acknowledge the very significant cultural role and 
reverence that the summit of Mauna Kea has al- 
ways had within the indigenous Hawaiian commu- 
nity. We are most fortunate to have the opportu- 
nity to conduct observations from this mountain. 
Radio data presented here were obtained using the 
Very Large Array of the National Radio Observa- 
tory, a facility of the National Science Foundation 
operated under cooperative agreement by Associ- 
ated Universities, Inc. 
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Table 2 
Extended Photometry Table of W40 IR Sources (Electronic Version Only) 



Index 



Q^J2000 



'.72000 



H 



K 



001 


277.762394 


-2.031111 


16.608 ± 0.145 


14.838 ± 0.203 


13.809 ± 0.272 


002 


277.764807 


-2.005333 


14.185 ± 0.144 


13.045 ± 0.203 


10.861 ± 0.272 


003 


277.767017 


-2.136254 


15.328 ± 0.144 


12.772 ± 0.203 


11.734 ± 0.272 


004 


277.767545 


-2.181531 


13.621 ± 0.144 


10.963 ± 0.203 


9.666 ± 0.272 


005 


277.767801 


-2.154991 


15.387 ± 0.144 


13.509 ± 0.204 


12.394 ± 0.272 


006 


277.767999 


-2.010032 


13.318 ± 0.144 


13.530 ± 0.203 


11.754 ± 0.272 


007 


277.769254 


-2.033301 


16.745 ± 0.146 


14.706 ± 0.203 


13.349 ± 0.272 


008 


277.770824 


-2.054937 


10.947 ± 0.144 


11.471 ± 0.203 


10.242 ± 0.272 


009 


277.770948 


-2.046558 


13.686 ± 0.144 


12.726 ± 0.203 


11.683 ± 0.272 


010 


277.771923 


-2.109518 


14.602 ± 0.145 


13.211 ± 0.203 


12.122 ± 0.273 


Oil 


277.772396 


-2.073557 


10.221 ± 0.144 


10.612 ± 0.203 


8.327 ± 0.272 



(complete table will be available in electronic version of AJ. Contact authors directly 
for access before publication.) 
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